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Abstract 



We study a new hierarchy of equations containing the Short Pulse equation, which describes the 
evolution of very short pulses in nonlinear media, and the Elastic Beam equation, which describes 
nonlinear transverse oscillations of elastic beams under tension. We show that the hierarchy of 
equations is integrable. We obtain the two compatible Hamiltonian structures. We construct an 
infinite series of both local and nonlocal conserved charges. A Lax description is presented for both 
systems. For the Elastic Beam equations we also obtain a nonstandard Lax representation. 



1 Introduction: 



The cubic nonlinear Schrodinger equation is used in the description of the propagation of pulses 
in nonlinear media such as optical fibers. Recently, technology progress for creating very short 
pulses was achieved, however, the description of the evolution of those pulses lies beyond the 
usual approximations leading to the nonlinear Schrodinger equation. Various approaches have been 
proposed to replace the nonlinear Schrodinger equation in these conditions. In Ref. Schafer and 
Wayne proposed an alternative model to approximate the evolution of very short pulses in nonlinear 
media. They derived the short pulse (SP) equation 

u xt = u+ l - {u 3 ) xx . (1) 

Chung et al 2 showed numerically that the SP equation provides a better approximation to the 
solution of Maxwell's equations than the nonlinear Schrodinger equation as the pulse length gets 
short. Also, Sakovich and Sakovich [2] have studied the integrability of |T]). 

In this paper we will study the integrability of the nonlocal representation of the SP equation 

© 

u t = («9 _1 n) + l;u 2 u x , (2) 

as well the hierarchy of equations associated with it. The equation in this form is more feasible for 
a Hamiltonian description and can be obtained integrating (^Q) with respect to x. In what follows 
we will refer to the equation (J2j) simply as the SP equation. 

The local nonlinear equation 

will also appear in our hierarchy of equations. It can be embedded in the Wadati-Konno-Ichikawa 
(WKI) system jl] and its x derivative can be shown to describe nonlinear transverse oscillation of 
elastic beams under tension [3]. Therefore, we will simply call © the equation for elastic beams 
(EB). 

Our paper is organized as follows. In Sec. 2, we shown that our system is integrable through its 
bi-Hamiltonian nature. We give the two Hamiltonian structures associated with the Hamiltonian 
description of equations (J2J) and ©. The method of prolongation is used to prove the Jacobi 
identity as well the compatibility of the Hamiltonian structures. In Sec. 3, we construct the 
recursion operator and its inverse to obtain the first local and nonlocal charges recursively. We 
also obtain the first local and nonlocal equations of the hierarchy recursively, which includes the 



(3) 
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EB and SP equations, respectively. In Se.c 4 we obtain the Lax representation for the system. For 
the EB equation, and for other local equations in the hierarchy, we also give a nonstandard Lax 
representation. In Sec. 5, we summarize our results and present our conclusions. 

The results in this paper are much like the ones obtained for the Harry Dym and Hunter-Saxton 1 
hierarchy of equations [B|, deformed Harry Dym and Hunter-Saxton hierarchy of equations [7j and 
for the Non-local Gas hierarchy of equations 8 . In these works we have as a main characteristic 
a hierarchy of integrable equations with positive and negative flows. Also, throughout this paper 
the calculations involving pseudo-differential operators were performed or checked by the computer 
algebra program PSEUDO [Sj. 



2 Bi-Hamiltonian Structure: 



Following 7, let us introduce 



F 2 ^l + ul, A^, (4) 



which satisfy the following useful properties 



F 2 (l - A 2 ) = 1 , 

J$ = AA x , 

^ = A x . (5) 



Using @ the SP equation © assumes the forms 

t) I -a 2 F 



2 

U 1 9 . 

At = F + 2 U Ax ' (6) 



and, similarly, the EB equation © can be written in one of the forms 

u t = A xx , 
Ft — AA XXX , 

At = ^. (7) 

1 In References |S| and |7j the equation obtained by Hunter and Saxton in SIAM J. Appl. Math. 51, 1498 (1991) 
was erroneously named, by the present author, Hunter-Zheng. We sincerely apologize the authors about this mistake. 
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Let us stress that the basic field is u and that A and F are just placeholders used to make expressions 
more compact. Also, it is interesting to observe that the EB equation when written in the form 
l [7| ) has the same structural form of the deformed Harry Dym equation studied in [7j, however, the 
definitions in Q are different. 

From (JBJ) and (JJJ) it is straightforward to note that 

H = -[dxF, (8) 



is conserved under both the SP and EB flows. 
Introducing the Clebsch potential 

u = <f>x , (9) 

the equation (J2J) can be written as 

<t>t = + \<Pl ■ (io) 

This equation (|l(Jj) can be obtained from a variational principle, 5 j dtdx C, with the Lagrangian 
density 

This is a first order Lagrangian density and, consequently, the Hamiltonian structure can be readily 
read out, or we can use, for example, Dirac's theory of constraints ^U] to obtain the Hamiltonian 
and the Hamiltonian operator associated with The Lagrangian is degenerate and the primary 
constraint is obtained to be 

<*> = 7T-^, (12) 

where tt = dC/d<f>t is the canonical momentum. The total Hamiltonian can be written as 
H T = f dx (rrcPt - C + A$) = j dx - + A U - ^ 



(13) 



where A is a Lagrange multiplier field. Using the canonical Poisson bracket relation 

{4>{x)^{y)} = 5{x-y) , (14) 
with all others vanishing, it follows that the requirement of the primary constraint to be stationary 
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under time evolution, 

{<S>(x),H T } = 0, 

determines the Lagrange multiplier field A in (|13ft and the system has no further constraints. 
Using the canonical Poisson bracket relations (|14j) . we can now calculate 

K(x, y) = {$(*), <%)} = ±d y S(y - x) - h x 5(x - y) . (15) 

This shows that the constraint (|12|) is second class and that the Dirac bracket between the basic 
variables has the form 

{4>{x),4>{y)} D = {4>{x)A{y)} - I dzdz'{<f>(x),<5>(z)}J(z,z'){<5>(z'),^y)} = J(x,y) , 
where J is the inverse of the Poisson bracket of the constraint (|15jl , 

dz K(x, z)J(z, y) = 5(x — y) . 



This last relation determines 



or 



where 



d x J(x,y) = -5(x - y) 



J(x,y) = V5{x - y) 



V = -d- 1 , (16) 

and can be thought of as the alternating step function in the coordinate space. We can now set 
the constraint (|12j) strongly to zero in (fl3|) to obtain 

H T = jdx - ^ . (17) 

Using © and the transformation properties of Hamiltonian operators |llj . we get 

v 1 = d (v) (dy = d , (is) 
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and the SP equation ((21) can be written in the Hamiltonian form as 

5H 2 



u t = T>i 



5u ' 



with H 2 given by 



H 2 



dx 



1 4 l/o-l \2 

24 U "2 (a U) 



(19) 



which can be easily checked to be conserved by both the SP and EB equations. 



We will show that the SP and EB equations belong to the same hierarchy of equations, at this 
point we note that 



u t = V] 

with Hq given by (JHJ), yields the EB equation j3J). 
It is easy to show that the charges 

1 



SH 
5u 



dxFAl 



1 



H\ = - I dxu 



(20) 
(21) 



are also conserved by both the SP and EB equations. Therefore, the SP equation (J2J) can be written 
in the Hamiltonian form as 



U t = V; 



SHi 

5u ' 



and the EB equation @ as 



ut = V; 



SH-i 
5u 



where we have defined 



V 2 = d 1 + u x d 1 u x = (F 2 - u x d 1 u xx ) d 



-1 



(22) 



The skew symmetry of this Hamiltonian structures is manifest. The proof of the Jacobi identity 
for this structure as well its compatibility with (|18j) can be shown through the standard method of 
prolongation ^2] which we describe briefly. 



6 



(23) 



We can construct the two bivectors associated with the two structures as 
Vl = \jdx {9AV 1 9} = ~fdx9A9 x , 

Q V2 = I J dx {6 A V 2 6] = ~ Jdx{9A (d^O) + u x 9 A (d~ V 9)} . 
Using the prolongation relations, 

pr v-Drfiu) = 9 X , 

prv V2 e(u) = (d~ 1 9) + u x (d~ 1 u x 9) , 
pr v V2 g(u x ) = (prv V20 (u)) x , 

it is straightforward to show that the prolongation of the bivector @t> 2 vanishes, 

prv V2 g(e V2 ) = , 

implying that T>2 satisfies Jacobi identity. Using ()23|) . it also follows that 

pvv Vl e (0© 2 ) + Pi"UD 2 e (OdJ = , 

showing that T>\ and T>2 are compatible. Namely, not only are T>\,T>2 genuine Hamiltonian struc- 
tures, any arbitrary linear combination of them is as well. As a result, the dynamical equations (J2J) 
and (pj are bi-Hamiltonian with the same compatible Hamiltonian structures and, consequently, 
are integrable ^3 EH • 



3 The Short Pulse Hierarchy: 

When a system is bi-Hamiltonian, we can naturally define a hierarchy of commuting flows through 
the relation 

Utn =K n [u]=V 1 6 -^=V 2 6 -^ , n = ... ,-2,-1,0,1,2,... . (24) 
du du 

For n = 1 and n = — 1 we get the SP and EB equations, respectively. For n = we have 

K = Vi^- = (du) = u x , 
ou 
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and 

K = V 2 ^- = ((F 2 - u x d- x u xx ) d-'Ax) 

= (F 2 A - Uxd^UxxA) = Fu x - u x F = , (25) 

where we have used ©, i.e., F 2 A = Fu x and u xx A = F x . Therefore, we would be lead to 
conclude that Hq is a Casimir of T>2- We can resolve this apparent contradiction being careful 
while performing calculations with the antiderivative <9 _1 . We can use the following representation 

/+oo 
dye(x-y)f(y), 
-oo 

where 

e{x -y) = 

Then, it can be shown that 



1/2 for x > y , 
T/2 for x < y . 



{d- l f x ) =f-\ (/(+oo) + /(-oo)) , (26) 

Now, if we assume — > as \x\ — > oo (which yields ^4|±oo = and -F|±oo = 1) then it follows 

(d- l A x ) = A and (d~ l F x ) = F — 1 . (27) 

Using this last result in the naive calculation l|25jl we obtain the desired term u x . This sort of 
missing or "ghost" terms given rise to apparent contradictions in nonlocal theories were already 
observed in the literature (see ^1] and references therein). 

Let us introduce the recursion operator following from the two Hamiltonian structures as 

R = V 2 V^ . (28) 

Then, it follows from (j2U) that 

SJ ^ = R^, n = 0,l,2,..., (29) 

where 

= d~ 2 + d~ x u x d~ x u x = d~ 2 (F 2 + u xx d' 1 ^) (30) 
is the adjoint of R. The conserved charges for the hierarchy can, of course, be determined recursively 
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from (|29j) and give the infinite set of (nonlocal) conserved Hamiltonians 



H n 



#1 



Ho 



dxF , 



- I dxu , 



dx 



dx 



J_ n 4_ \g-i u) 2 
24 2 y ' 



720 2 V ' 6 V ; 4 V 7 



; . (si) 

The corresponding flows (the first few, since the equations become extremely nonlocal as we proceed 
further in the recursion) have the forms 

Uf = u x , 



u tl = (d 1 u) + -u 2 u x , 

u t2 = (d~ 3 u) + i(3-V) + u x (d- 1 (u x {d~ 2 u))) + i 



(32) 



24 



U Ur, 



For negative values of n, the gradients of the Hamiltonians will satisfy the recursion |3] 

» — 

Writing the recursion operator (|28[) in the form 

i? = (f 2 - n x a- 1 ^) a- 2 , 

and using the identities © the inverse can be easily checked to be 

I / ' , /I /t 1 _fl2 1 Q 77/5-1 Uxx 



R-' = d<[- 2 +Ad->A x =d<— dFd F3 



(33) 



(34) 



(35) 



Note that (|35|) can be recognized as the recursion operator obtained in [3] for the SP equation 
in the form (^Q) using cyclic basis techniques. The corresponding conserved charges can now be 
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recursively constructed from (|33j) and have the forms 
H-i = \j dxFA 2 x , 

1 f / A A 3 A 2 A 2 A 2 \ 

H - 3 = i6 J dx [ FA * + 8_ f ££ - 12 + ) ' 

i • (36) 



The corresponding flows (the first few) have the forms 
ut_i — A xx , 

(37) 



_ 1 A xxxx 1 A 2 A XX ^ A XXX A X A 3 A 2 X A 2 A 2 , 3 ^ 1 

u '3-l f 4 + 2 F 2 F 2 2 F 2 ^Ac^ 1 + o A x A 



XX 



4 The Lax Representation: 

Conserved charges for our systems can be determined in principle recursively from 1)29(1 and (|33|) . 
However, to construct the conserved charges directly we look for a Lax representation for the system 
of SP and EB equations. 

It is well known ^3 EH that for a bi-Hamiltonian system of evolution equations, ut n = K n [u], 
a natural Lax description 

is easily obtained where, we can identify 

M = R, 
B = K' n . 



10 



Here R is the recursion operator (|34|) and K' n represents the Frechet derivative of K n , defined by 

K' n [u] v = 4" K n [u + ev] 

de e=o 

For the SP and EB system of equations in @ and @ respectively, we have 

B SP = K'^d-' + ^du 2 , 

B EB = K'_ 1 =d 2 F~ 3 d. 



The two systems have the same M = R given in (|28|). It can now be checked that 

dM 



at 



[B SP ,M] , 



™ = [B™,M], (38) 

do indeed generate the SP and the EB equations and, thereby, provide a Lax pair for the system. 

A Lax representation directly gives the conserved charges of the system. From the structure of 
Q38JI . it follows that TrM^^ are conserved, where "Tr" represents the Adler's trace |17j . We note 
that 

2n+l 

TrM^ - = , n > 1 
Trikfi = dxF , 



TrM~5 = ^ J dxFAl , 
TrM-l = - ! dx ( FAi - 4 ^ 



(39) 



These charges correspond (up to multiplicative constants) to the ones given in (|3"H|) . constructed 

_ 2n— 1 

earlier by recursion. In fact, all H^ n with positive n > can be constructed from TrM 2 and 
by construction (namely, because of the nature of (|38|l). they are conserved under both the SP 
and EB flows. However, as is clear from ()39|) . this procedure does not yield the nonlocal charges 
H n with positive integer values. The construction of these nonlocal charges relies primarily on the 
recursion relation ([29)1 . As in our papers [01 HI IE] it remains an interesting question to construct 
these charges in a more direct manner. 

We can also obtain a Lax representation, in a Gelfand-Dikii form, for the EB hierarchy of 
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equations, in a straightforward way. It is well known (see [7| and references therein) that the 
recursion operator defines a eigenvalue problem of the type 



(1 - \R}) <j) = 



(40) 



for the eigenfunction with eigenvalue 1/A. Since we know the inverse of the recursion operator 
we obtain from ()4(Jj) 



= 0, 



(41) 



which defines a Lax eigenvalue problem. Using (|35|) we identify the Lax operator for our system to 
be 

R ] Y=jdjd + A x d- l A x d. 

It can be readily checked that the hierarchy of EB equations (|37|) (up to multiplicative constants) 
can be obtained from the non-standard Lax equation 



dL 

dt- n 



( L (*»+l)/2) £ 



n = 1,2,3,. 



The conserved quantities (jSHJ) an d Hq in IfMljl for this system (up to multiplicative constants) can 
be obtained from Tr L^ 2 ™" 1 )/ 2 , n = 0, 1,2, . . .. 



5 Conclusion: 



In this paper, the Short Pulse and Elastic Beam equations were shown to belong to the same 
hierarchy corresponding to positive and negative flows. We have shown that these system are bi- 
Hamiltonian and using recursion we construct infinite series of both local and nonlocal conserved 
charges as well the respective hierarchy of equations. A Lax pair for the Short Pulse and Elastic 
Beam equations was derived. The Lax operator yielded the local charges via Adler's trace, however, 
the construction of the nonlocal charges through this Lax operator is unknown by us. Also, a 
nonstandard Lax representation in a Gelfand-Dikii form for the Elastic Beam hierarchy of equation 
was described. 
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